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Bond-disordered Anderson model in two dimensions on a square lattice 
is studied numerically near the band center by calculating density of states 
(DoS), multifractal properties of eigenstates and the localization length. DoS 
divergence at the band center is studied and compared with Gade's result 
[Nucl. Phys. B 398, 499 (1993)] and the powerlaw. Although Gade's form 
describes accurately DoS of finite size systems near the band-center, it fails 
to describe the calculated part of DoS of the infinite system, and a new ex- 
pression is proposed. Study of the level spacing distributions reveals that the 
state closest to the band center and the next one have different level spacing 
distribution than the pairs of states away from the band center. Multifrac- 
tal properties of finite systems furthermore show that scaling of eigenstates 
changes discontinuously near the band center. This unusual behavior suggests 
the existence of a new divergent length scale, whose existence is explained as 
the finite size manifestation of the band center critical point of the infinite 
system, and the critical exponent of the correlation length is calculated by a 
finite size scaling. Furthermore, study of scaling of Lyapunov exponents of 
transfer matrices of long stripes indicates that for a long stripe of any width 
there is an energy region around band center within which the Lyapunov ex- 
ponents cannot be described by one-parameter scaling. This region, however, 
vanishes in the limit of the infinite square lattice when one-parameter scaling 
is restored, and the scaling exponent calculated, in agreement with the result 
of the finite size scaling analysis. 

PACS numbers: 71.23.An, 72.15.Rn, 73.20.Fz 



I. INTRODUCTION 

A quantum particle moving in a random potential undergoes Anderson localization quan- 
tum phase transition in three dimensions with increasing of the strength of disorder |T]- [J . 
The order parameter characterizing the localized phase is the inverse localization length 
||, describing the exponential decay of the envelope of eigenstates. When the criti- 
cal point is being approached from the localized phase, localization length, which depends 
for a given energy only on the strength of the disorder, increases with decreasing of the 
disorder strength and finally diverges as a power-law at a particular disorder strength. Fur- 
ther decrease of disorder strength then makes the eigenstate extended throughout the whole 
system. Simultaneously, on length scales smaller than the localization length, eigenstates 
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exhibit multifractal scaling behavior characterized by anomalous scaling of the inverse par- 
ticipation numbers (for definitions and references, see Sec. 0). 

This basic phenomena, together with the work of Thouless on the scaling of conductance 
in finite-size systems ||, led to the scaling theory of localization one of the main con- 
sequences of which is the absence of extended states in two dimensional disordered systems, 
with two dimensions being the lower critical dimension of the transition. If the spin-orbit in- 
teraction is present, however, picture changes and two dimensional systems from symplectic 
ensemble exhibit localization transition even in two dimensions as opposed to systems from 
orthogonal ensemble which have all states localized . Presence of the strong magnetic field 
in two dimensional disordered systems, on the other hand, leads to a completely different 
behavior - the integer quantum Hall effect (IQHE) - where critical states are present at the 
middle of each of disorder-broadened Landau levels ||. 

Another class of models exhibiting localization properties different from the systems 
mentioned above are systems with chiral (particle-hole) symmetry. Such systems are defined 
on a bipartite lattice with only hopping (off-diagonal or bond-) disorder. Wegner was first 
to realize the importance of this symmetry in disordered systems [|H3|-|T5|j, and even one- 



dimensional systems with this symmetry are known to have peculiar properties, such as 
diverging DoS at the band center |TD[, where the eigenstate decays as exp(— jy/r) |TT| . |T2 



in contrast to one- dimensional site-disordered systems which have DoS bounded [16| and all 
states localized. 

There are several models with chiral symmetry that have been extensively studied. The 
simplest two, in the sense that only one orbital per site and nearest neighbor hopping are 
included, time-reversal symmetry present and the spin not relevant, are the Anderson bond- 
disordered model (ABD) [[1C 12, 1^,20] and the random Dirac fermion model (RDF) |17,18 



The main difference between these two models is that, in the non-disordered case, ABD 
model has a line of points as the Fermi surface at half-filling while RDF has a point Fermi 
surface and linear dispersion of energies. 

This work is concerned with the ABD model on a square lattice of size L and periodic 
boundary conditions, defined by the Hamiltonian: 



H = ~ e o E {hi c h + H - c ) > (!) 

where brackets denote neighboring sites on the lattice, q is annihilation operator of the 
electron at site i, and £'s are uniformly distributed random variables tij G (1 — 2w, 1), with 
< w < 1. They represent random hopping energies between nearest neighbors, expressed 
in units of energy €q, which is set to 1 hereafter. 

Interest in this model mainly comes from its unusual scaling properties at the band center, 
where Soukoulis et. al [HJ have found critical state at the band center using Green's function 
H and transfer matrix method (TMM) j2!|. More recent TMM calculation by Eilmes et. al 
20| confirmed this result with a higher accuracy and showed the validity of one-parameter 



scaling not too close to the band center. Nevertheless, Miler and Wang have found in 
their study of two models with chiral symmetry an apparent band of extended states near 
the band center, and it remained unclear what is the fate of these states in the infinite 2D 
system. Another study |2B| yet showed that scaling exponent of the average participation 
number changed discontinuously near the band center, and the explicit dependence of this 
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energy on the system size proposed by authors implied the existence of another diverging 
length scale in the problem. The last effect is rather subtle to calculate and led to a different 



participation number scaling exponent of the ABD model at the band center in Ref. GO 



compared to the one calculated here, as discussed in detail in Sec. [V] below. Furthermore, 
Brouwer et al. |23| have calculated conductance distribution of quantum wires described by 
(HI), and showed its non- universality and necessity to introduce an additional microscopic 
parameter. 

It is thus goal of this paper to present a detailed study of the scaling of localization 
length on the approach to the band center for an infinite 2D square lattice, and test the 
validity of one-parameter scaling, as well as to calculate the multifractal properties of the 
electron probability density on length scales smaller than the localization length. Also a new 
analytical expression for the DoS of the infinite two-dimensional system near the band center 
is proposed. This paper is organized as follows: Some general properties and exact results are 
presented in Sec. |T[ Calculation of DoS is presented and analyzed in Sec. [TTI|; Sec. [TV] contains 
analysis of level spacing distributions between the nearest neighbors; Multifractal properties 
of eigenstates are studied in Sec. |V|; Scaling of the Lyapunov exponents of transfer matrices 
of long strips and the scaling of localization length are studied in Sec. |VI|; and, finally, Sec. 
V 111 summarizes results of this work. 



II. SOME GENERAL PROPERTIES OF LATTICE HAMILTONIANS WITH 

CHIRAL SYMMETRY 

Suppose that the lattice is composed of two sublattices A and B with, respectively, Na 
and Nb sites. The corresponding bond-disordered Hamiltonian with chiral symmetry then 
has the form: 

H= E M^ + Hx), (2) 

ieAjeB 

It is easy to show that for every eigenstate with energy E there is an eigenstate with 
energy —E with a wavefunction that has the opposite sign at each site of one of the two 
sublattices. 

If the total number of cites N = Na + Nb is odd and open boundaries condition is 
applied (in order to keep the symmetry), then, since all eigenstates come in the opposite 
energy pairs, there will be exactly one state with eigenenergy 0. This can be furthermore 
generalized, and if m = Na — Nb > 0, there exist exactly m zero-energy eigenstates that 



have vanishing amplitude on the sublattice B fl5||2 



On the other hand, if m = 0, the electron has equal probability of occupying each of 
the two sublattices. To show this, (0) is represented in the basis where the first and second 
half of basis vectors are eigenstates of the position operator on sites of sublattice A and B, 
respectively. The Hamiltonian is then represented as 

"(if). (3 » 

where M is a square matrix of hopping elements from one sublattice to the other. Eigenstate 

L/A _/W 



Mb) 



satisfies: 
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E={^\H\^) = 2Re(^|M|^ B ). 



(4) 



On the other hand, 

From @) and (H) now follows that, for ABD model, (i^aI^a) — 1/2- 

In this work only even L finite size systems on a square lattice with periodic boundary 
conditions are studied, because one of the main goals of this work is to understand the 
vicinity of the critical point of ABD model on the infinite square lattice, which in turn has 
m = 0, while the limit L — > oo for odd L and open boundary conditions has m = 1. 



III. DENSITY OF STATES NEAR THE BAND CENTER 

Density of states (DoS) is calculated by exact numerical diagonalization of finite size 
Hamiltonians for various L for many configurations of disorder, and binning of eigenenergies. 
Obtained DoS for each system size, Pl(E), are normalized to 1. The L dependent parts of 
such obtained Pl{E) are then removed leaving L independent DoS p(E), which is therefore 
expected to be correct in the L — > oo limit. The removal of finite size dependency is based 
on an observation that DoS converges quickly away from the band center with increasing 
of L. Thus, only a small number of eigenenergies (up to 20) closest to the band center and 
corresponding DoS histograms has been calculated for each L. The calculated Pl(E) plotted 
on a single graph revealed that three bins closest to the band center are where the system 
size dependence sets in. Their removal thus led to the DoS p{E) of the infinite system. 

Results for pi(\E\) are given in Fig. [I], for system sizes L = 10, 20, ...,60 and number of 
disordered configurations ranging, respectively, from 160000 to 4100. They are fitted to a 
powerlaw divergence Pl(E) = Ci,\E\~ aL , as well as to the Gade's result fl4fl : 



p L (E) = C L j^ exp (-K Ly f-ln\E\) . (6) 

All the calculations were done for several different number of bins, and the obtained values 
of fitting parameters were the same within error bars. 

The figure [l] shows that (|6]) describes p^ very accurately for L > 40, including the size- 
dependent part. The powerlaw, on the other hand, also describes the data accurately for 
same system sizes, but fails to describe the L-dependent part of p^. Despite this, neither of 
the two forms describe the whole p(E) accurately. Instead, the expression found to best fit 
the obtained L-independent DoS, given in Fig. ||, is 

p(E) = C^- exp (-«y-ln|£|) , (7) 



with k = 1.345 ± 0.005 and C = 1.30 ± 0.03, represented by the full line in the same figure. 
The observed range in which (0) is accurate is for all the energies studied smaller than 
6 x 10" 2 . 
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IV. DISTRIBUTION OF THE NEAREST-NEIGHBOR LEVEL SPACINGS 



In the localized regime, an eigenstate is determined mainly by a local configuration of 
disorder where the wavefunction is localized, and two eigenstates close in energy are spatially 
far apart. Level repulsion is therefore absent and the distribution of the nearest neighbor 
level spacings s = E i+ i — E { is Poissonian [EH] , 



D P (s) = ^exp(-£). (8) 

where 6 = (s) is the mean level spacing. 

In the delocalized phase, on the other hand, eigenstates are extended throughout the 
system and level repulsion becomes significant for eigenstates with close energies. In the 
infinite 3D Anderson site-disordered model (ASD), 

H = e i n i ~ { c \cj + , (9) 

with uniformly distributed e.; G (— W/2, W/2), distribution of level spacings becomes that of 
Gaussian orthogonal ensemble (GOE) [pHJ, very accurately described by the Wigner surmise 



JM.) = ^«p (-£(£)). do) 

In finite-size systems, localized states are on average at a distance L rather than infinitely 
far apart. This leads to a repulsion between adjacent energy levels and non- universal dis- 
tribution Dl(s) . Shklovskii et. al. have shown that -Dl(s) of the 3D site-disordered 
Anderson model exhibits linear dependence on s characteristic for Dyy(s) for small s and 
exponential tail characteristic for Dp(s) for large s. They were able, from the finite size 
scaling analysis of the tail, to accurately determine the critical point and exponent. In the 
infinite size limit, they have recovered not only Dp(s) in the insulating phase and D\y(s) 
in the conducting phase, but also a system size independent non-universal distribution at 



the critical point which was furthermore shown by Braun et. al. [37[] to be dependent on 
boundary conditions. This method was also used for an accurate determination of localiza- 
tion length in two dimensional ASD model P%[|, confirming the absence of delocalized states 



following the scenario of the insulating phase from Ref. |26| described above. 

To see the effect of symmetry of the Hamiltonian (^) on the distribution of level spacings, 
let us for a moment consider i-th eigenenergy Ei of the ASD model. Upon averaging over 
disorder, the Ei will be distributed between E™ n and E™ ax according to some distribution. 
Some of the eigenenergies, for i close to N/2, will have E™™ < < E™ ax . This is, however, 
forbidden for eigenstates of the ABD model since every Ei of @) is negative for i < N/2 and 
positive for i > N/2. This means that eigenenergies of (|2]) close to the band center are effec- 
tively pushed away from it due to the symmetry. If L is much smaller than the localization 
length, states will be repelled among themselves due to their large spatial overlap. But the 
two states closest to the band center, being simply related to each another by the symmetry, 
will not repel at all, i.e. the state closest to the band center is at the (high energy) end of 
the spectrum. Thus, these two states are distributed around zero, where distributions of all 
other individual levels go to zero. 
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To explore consequences of this simple analysis, the level spacing distribution is calcu- 
lated between each pair of adjacent levels separately. Let us denote by D^s) level spacing 
distribution between % and i + 1 energy level after unfolding of spectrum P5[ , i.e. expressing 
level energies in units of mean level spacing, where < E\ < E 2 < ... < E N / 2 . Figure || 
shows -Di(s), D$(s), for L = 20, 40 and, respectively 150000, 120000 configurations, and it 
can be seen that -Di(s) is distinctly different than D 2 (s), Z?s(s). The same effect was also 
present for L = 10, 30, while for L = 50 and 60 the number of disorder configurations was 
insufficient for accurate enough determination of individual D^s) ||30|| . This illustrates how 
the presence of chiral symmetry can profoundly influence spectral characteristics near the 
band center, despite the fact that DoS of ABD and ASD models seem to have same shape 
for adequately chosen pairs of disorder parameters w and W away from the band center 
(and after rescaling of 6q) JHJ. 



V. MULTIFRACTALITY OF EIGENSTATES 

Eigenstate of an electron in the random potential fluctuates from site to site and it was 
proposed that the eigenstate at the mobility edge in disordered systems in general should 
have fractal structure pi , and shown that even localized states in one and two dimensions 



exhibit fractal character on length scales smaller than the localization length [|32"| , |3"3" |. 

Inverse participation numbers Z q (IPN) are particularly convenient quantities to describe 
scaling properties of probability distribution of the electron. IPN of an eigenstate \l/ are 
defined as: 

Z g (*)= Sl^rOI 3 ". (11) 

i=l 

Intuitively, their meaning can be seen by looking at the participation number Z 2 (^)~ 1 : it 
is equal to 1 for a state localized at one site and to N for plane-waves. Participation number 
thus gives generally the number of sites at which the wavefunction is significantly different 
than zero. Participation numbers Zq^)^ 1 then generalize this by giving the number of sites 
where probability distribution of electron is very high (for large positive g's), very low (for 
large negative g's), and in between these extrema, continuously parameterized by q. 

More convenient, with the advantage of being defined as averages over disorder at a given 
energy E, are IPN defined as functions of E and system size L, 

Z q {L,E) = (Z q ^)8{E^)-E)), (12) 

where the brackets denote averaging over disorder. Z q (E,L) can be numerically calculated 
by averaging ( |TT| ) over all eigenstates from M configurations of disorder belonging to an 
energy interval of width AE around E, and studying the limit AE — > for large M [^3] . 



Wegner pioneered this kind of investigations, and Castellani and Peliti proposed 
that eigenstates near the critical point are multifractal on length scales smaller than £. The 
most important feature of IPN of eigenstates is their scaling with system size and energy 



34,35 : 



Z q (L,E)~L- T \ (13) 
Z 9 (L,Er)~\E-E c \**, (14) 



6 



where E c is the critical energy. The former scaling is present at any E for L <C £,{E), while 
the latter holds in the critical region of the transition . 

Within the framework of multifractality [^,^], electron probability density is charac- 
terized by several quantities that can be derived from r q — the generalized dimension D q 
and the singularity strength a q of the g-th singularity with the fractal dimension f q : 

dr 

(q - l)D q = r q , a q = -j± , f q (a q ) = a q q - r q . (15) 

D q represents generalization of the fractal dimension, and it is constant and equal to the 
fractal dimension for ordinary fractals, while f q {a q ) is the singularity strength spectrum 
describing multifractal as an interlaced set of fractals with fractal dimensions f q , where the 
measure on the g-th fractal scales as a powerlaw with exponent a q . These quantities have 
several general properties: D is the fractal dimension of the support (two in this work); Di 
is called information dimension since it describes scaling of the entropy of the measure J38| . 
and there exist finite D min = and D max = _D (? ^_ 00 . 

This work is concerned mainly with the spectrum of generalized dimensions D q charac- 
terizing the spatial structure of eigenstates on the length scales smaller than the localization 
length, while the properties of ir q will be discussed elsewhere. Before detailed discussion of 
the results, an overview of the main results of this part of the paper is given. Scaling of 
IPN at the band center is calculated first, and shown that scaling properties (that is, whole 
spectrum of generalized dimensions D q ) changes discontinuously near the band center, at 
an energy E'(L) for a range of L studied. It is then shown that E' can be quite accurately 
identified with the half of the width of the energy range around band center within which 
two states occur on average in the ensemble of disordered systems. Existence of this energy 
reveals the existence of a length scale diverging when E — > 0, that is the system 

size at which IPN change their scaling properties from one powerlaw dependence on L to 
another. This change is then explained as a finite size manifestation of the critical point, 
and the critical exponent calculated by a finite-size analysis. 

Calculation of Z q (E, L) starts with calculation of Z q (E, L, AE), which is just IPN aver- 
aged over all eigenstates from an energy interval {E — AE/2,E + AE/2), taken from ]Vq 
realizations of disorder, followed by studying the limit AE — > ||23|| . Results at the band 
center for system L = 80 and several different g's, are presented in Fig. |j. The error bars in 
the figure are taken to be the standard deviation of average value. 

The figure suggests the existence of an energy E' independent of q (and therefore defined 
by the whole multifractal measure) such that decreasing AE below E' does not change 
Z q (E, L, AE) significantly. Decreasing of Z q to a smaller extend, however, is still present for 
AE < E', and the main source of this is the mismatch between average and typical value 
of IPN at a given energy. Thus, the effect should become smaller as the number of disorder 
configurations that are averaged over is increased, and Z q (E, L, AE < E') ~ Z q (E, L) up to 
the corresponding statistical error. This can be seen in Fig. [| and |5], where values of Z q for 
AE < E' are approximately constant within the error bars (as indicated by the horizontal 
dashed lines), while for AE > E' there is approximately linear dependence of Z q on the bin 
size AE. In this sense Figure [5] suggests that E' exists for all the systems studied (and can 
be shown to be independent of q for each of them analogously as shown in Fig. |5|). 

Analogous analysis is carried out for energies away from the band center and Z q (E, L) 
determined accordingly, where it turns out that the convergence for these energies is slightly 
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easier to establish and occurs at larger AE than at the band center. Such obtained energy 
intervals used in calculations of IPN for different E were small enough so that there was 
practically no overlap among them. 

The vertical dashed lines in Fig. (| and Fig. |5] represent half of the energy E 2 (L), defined 
as the width of the energy interval around zero in which every system from the ensemble of 



disordered systems has two states on average [23 
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1 = L 2 p L (e)de. (16) 

Je 2 /2 

From the figures one can see that E' ps E 2 /2 for all system sizes studied except L = 40 (the 
smallest system studied, not shown in Fig. (5|), where the convergence seems to be somewhat 
slower. Equation ( |T6"D defines a new length scale (,'(E) that can be described as the system 
size L for a given energy E such that the number of states within the energy interval between 
— E and E is two on average. This can be defined as 

e(E) = E 2 \E), (17) 

where E 2 1 (E) is the inverse function of the E 2 (L), which is defined by fli~6"l) . Since E 2 (L) 
goes to zero when L — > 00, the new length scale diverges when E — > 0. 

It is tempting to integrate results for p L from Sec. |T| to obtain £'(E) explicitly. This, 
however, does not give the correct result since the whole analysis of DoS from Sec. |TTTj is 
done for energies larger than the width of the distribution of two states closest to the band 
center, which in turn defines £' in (|T7|). In other words, there is an energy cutoff, vanishing 
when L — > 00, below which the fits are not accurate, most obviously seen by noticing that 
all of the assumed analytical forms of p^ are diverging at the band center, while the actual 
Pl is not. 

Nontrivial feature connected with the existence of the new length scale £' is that scaling 
exponents r q {E) are different for L < £'{E) and L > £'(E). This is a generalization of 



findings from Ref. [p3|| , where a deviation from the powerlaw scaling of the average par- 
ticipation number for different AE at the band center appeared whenever AE exceeded 
E'(L) in several models with chiral symmetry. It can be straightforwardly shown that same 
happens with scaling of Z q (E = 0,L,AE). This suggests that a new scaling characteristic 
should be attributed to the two states closest to the band center, and that E' scales as E 2 , 
with a coefficient of proportionality close to one. This also implies that the corrections to 
the constant Z q (E, L, AE) for AE < E' are small. Therefore, the length at which change 
of scaling occurs should be close to and depend on the energy proportionally to £'(E), while 
the change of scaling should be a narrow crossover, as opposed to much broader crossover 



from powerlaw to constant IPN that occurs at the termination of multifractal scaling ([13 
for L ps £. 

This is compatible with the results presented in Fig. || which gives calculated Z q (L, E), 
for q = 0.9, 2 as well as for all L and E studied. Smaller g's allow for more accurate 
determination of the scaling properties, and it can be seen from the data that for both q's 
there are two characteristic scaling behaviors — one occurs at the band center and nearby 
energies for smaller L, while the other scaling holds at energies further away from the band 
center, and for larger L at energies close to the band center. 
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Scaling exponents r q are determined form the linear regression of the data for system 
sizes L = 50,60, 100. Goodness of the powerlaw fit is quantitatively characterized by a 
coefficient j q (E) next to the quadratic term from an additional quadratic fit of the data. 
Such obtained r q and 7 g from the data in Fig. || are presented in Fig. [F[ Results show the 
existence of three different cases: (i) away from the band center, 7 g « and r q is independent 
of E, indicating powerlaw dependence of IPN on L; (ii) approaching the band center, 7 g 
becomes different than zero, indicating that powerlaw is not obeyed. This is due to the 
emergence of the new scaling for system sizes L > £'(E) discussed above and present in 
Fig. and (iii) for E = powerlaw is obeyed again (7^ « 0), but with a different r q than 
for energies away from the band center. 

Discontinuity of r q (E) for all the other q's studied naturally leads to two different spec- 
tra of generalized dimensions, and Fig. |8| shows calculated D q for all energies except the 
three nonzero energies closest to the band center (which cannot give r q from the fitting 
procedure used here due to the change of scaling properties discussed above). In particular, 
the participation number grows with the number of sites L 2 as a powerlaw with exponents 
(3(E = 0) = 0.25 ± 0.02, and (3(E ^ 0) = 0.55 ± 0.05. This should be compared with 
the result for ABD model of Ref. 0, @(E = 0) = 0.50 ± 0.06. It is easy to explain this 
discrepancy since the bin size AE = 4 x 10~ 4 used in Ref. [2~0"| was, depending on the system 



size, roughly over an order of magnitude too large to detect the correct scaling behavior, 
and therefore /3(E 7^ 0) was obtained instead. 

D q is calculated for only three negative values of q, mainly because IPN for negative g's 
are determined mostly by the parts of eigenstates with the smallest probability to find the 
electron, which in turn acquire the highest relative error during numerical diagonalization 
of the Hamiltonian. Difficulties in calculating D q in this regime even arose suspicion that 
multifractality might break down for negative q [ID]. It is thus important to show that D q 



is defined for negative g's as well as for positive ones. Accuracy of all the calculations of 
IPN done in this section can be straightforwardly improved by increasing the number of 
configurations of disorder that was averaged over. This would lead to smaller error bars 
of all the quantities calculated, as well as to the wider range of g's for which D q can be 
calculated. 

Results of this section give the following picture of the scaling of IPN with system size: 
for any energy E close enough to the band center, there exist two powerlaw scalings: one 
for L < L < (,'(E) described by the set of exponents T q (E = 0), and another one for 
£'(E) ^5 L <C £(E), described by a different set of exponents r q {E 7^ 0), which leads to the 
two different spectra of generalized dimensions D q . 

Dependence of the new length scale on energy, £'(E), can be easily determined from QT7D 
by integrating the actual numerical data for pl(E), and the result, obtained from Fig. |9|, 
gives £,'{E) oc \E\~ V ', with 

v = 0.35 ± 0.01 . (18) 

The meaning of this new length scale and corresponding exponent v' can be understood 
by assuming that the additional scaling of the two states is due to the finite size effect 
of "smearing" of the E = critical point of the infinite system, because exactly the two 
states closest to the band center become critical when L — > 00. The critical energy then 
changes by Ae which is, in a finite system of size L (in units of the lattice spacing), equal 
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to AEziL) oc L~ 5 = L~ x l v ' . If furthermore £ oc \E\~ y in the infinite system, the shift Ae c 
of the critical energy in the system of size L is Ae c oc L -1 ^, from the general theory of 
the finite size scaling PH . Therefore, i/' equals u, the critical exponent of the correlation 
length of the infinite two-dimensional system. It should be noticed that the finite size scaling 
applied here is somewhat different than usual, where Ae c = \E — E c \/E c Here, E c = 
and Ae c = \E — E c \, where both energies are expressed in units of e , as discussed in the 
Introduction. Thus, energy eo appears in the denominator of Ae c rather than the critical 
energy itself. 



VI. LOCALIZATION LENGTH NEAR THE BAND CENTER 

In order to calculate the localization length of the infinite two dimensional system, finite- 
size scaling analysis of MacKinnon and Kramer M (FSS) is applied to TMM of Pichard and 
Sarma ||22|| . In this analysis, inverse Lyapunov exponents (ILE) of the transfer matrix of the 



ABD model [^0] : ^4| of a long quasi one dimensional strip of width M are calculated for several 
energies near the band center and one parameter scaling analysis applied to the largest ILE, 
from which the correlation length of 2D system is calculated. The scaling analysis consists 
of assuming that the change of the largest ILE, A(E, M), due to rescaling M — > bM can be 
compensated by an appropriate change of energy, after which A will remain the same, which 
implies that || 

A(E,M) = A(M/£(E)). (19) 

Figure [U] shows the calculated A(E, M) for several energies E close to the band center and 
for various strips up to 128 sites wide. The figure also gives the second largest (dashed line) 
renormalized ILE, A 2 , for the two lowest non-zero energies studied (E = 10~ 5 and 10 -6 ) and 
for M < 20. All values are obtained with relative error of 1% or better. 

At E — and for M even, all ILE become doubly degenerate due to the presence of 



chiral symmetry pl| , p4| , and they scale linearly with M for M > 16, reflecting the scale 
invariance of A characteristic of a critical state. To see the effect this degeneracy has 
on scaling properties of ILE, we should recall that ILE of transfer matrices of disordered 



systems repel each other in general []42| and become self-averaging quantities for sufficiently 
long stripes |fl2| . The symmetry now enforces "dimerization" of pairs of ILE, acting as an 
effective attractive force between each pair of ILE that become degenerate at zero energy, 



as in Fig. ITU], where, for a fixed M, A and A 2 are closer together for the smaller energy. 



The largest ILE thus decreases in a strip of width M on approaching the band center, as 
in models with chiral symmetry studied in Ref. [|2TJ. On the other hand, at any given 
energy, every pair of ILE becomes more repelling with increasing of M. Such increased 
number of ILE thus diminishes the attractive effect of the symmetry, and, depending on the 
relative strengths of the two effects, there are two regimes: (i) for smaller M, the largest 
ILE increases approximately linearly with M, and, since the slope of the rise changes with 
energy, FSS cannot be done; (ii) for sufficiently large M, on the other hand, the repulsion 
among ILE due to disorder dominates and FSS is possible. 

The one-parameter universal function A(M / £(E)) , obtained for system sizes M > 50 
for energies \E\ > 10~ 4 , and for system sizes M > 64 for energies 10~ 5 < \E\ < 10~ 4 , is 
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presented in Fig. |TT[ The obtained localization length £(E), in units of the localization 
length of the smallest energy studied (E = 0.1), is shown in the inset of the same figure. 
Calculated £(E)/£(0.1) is fitted to the powerlaw for energies 10~ 5 < E < 10 -3 , and the 
result suggests powerlaw diverging localization length at the band center, with the exponent 

v = 0.335 ±0.034, (20) 

in agreement with the result ( |18D of Sec. [V|. 

In Ref . [GDI , the authors also calculated one-parameter scaling curves for energies further 



away from the band center than in this work and found that there is a universal curve for a 
large set of ABD systems with different w that they studied. The obtained one-parameter 
function suggested that A(x) either grows indefinitely when E — > 0, which is similar to the 
two dimensional ASD model || or becomes constant. If this result were correct, however, for 
the smaller energies than studied there, it would imply that localization length at the band 
center either diverges logarithmically, as in the 2D ASD model, or that it has a finite limiting 
value. We see therefore that, although correct in the energy range studied in Ref. [j2U[ |, the 
picture of a single scaling curve for different w from the same work should break down for 
smaller energies, since the powerlaw diverging £ requires that din A/ din M goes to zero 
when E — > 0, and separate one-parameter scaling curves exist for each w. 

Some additional analysis of the results can be done by introducing A max (M), defined as 
the maximal A(E, M) for a given strip width M. Importance of this quantity comes from 
the fact that points where A reaches its maximum for various energies cannot be described 
by FSS, and, at the same time, A max limits by its definition possible values that A(x) 
(obtained from the FSS) can have. The main observation is that A max seems to grow slower 
than linearly in Fig. [II]. Linear (or faster) growth would imply existence of an additional 
energy scale below which FSS would break down for all large M and A(E, M) would grow 
linearly and indefinitely with M in the same figure, implying the existence of a whole band 
of extended states where one-parameter scaling would not hold. Slower than linear growth 
of A max seen, therefore, furthermore suggests the existence of a single critical point at the 
band center of the system on an infinite square lattice. 



VII. CONCLUSIONS 

In conclusions, main results of this work are summarized. Density of states of finite 
size systems are calculated and shown the validity of Gade's expression (|6]). The calculated 
part of the DoS for the system on an infinite square lattice, however, suggests a different 
dependence on energy near the band center, 

p(E) = C-^- exp (-^-ln\E\) , (21) 

with k = 1.345 ± 0.005 and C = 1.30 ± 0.03. 

Other calculated quantities share in common a qualitative feature of a discontinuous 
change near or at the band center. The nearest neighbor level spacing distributions, for 
instance, between the state closest to the band center and the next one seem to be distinctly 
different than the level spacing distributions between other neighboring states. This was 
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argued to be connected to the chiral symmetry of the model, which places the two states 
closest to E = at the (high energy) end of the spectrum. These two states are furthermore 
found to play a crucial role in explaining the discontinuous change of the scaling properties 
of IPN near the band center. Extrapolation to the limit of infinite square lattice then led 
to the two different spectra of generalized dimensions D q — one for E ^ (present at the 
length scales smaller than the localization length), and another one for E = (present at all 
length scales). Finite size scaling associated with the effects that finite L has on the critical 
band center states of the infinite system led to the value v = 0.35 ± 0.01 of the critical 
exponent of the localization length. 

Additional scaling analysis of ILE of transfer matrices of long quasi- ID systems gave the 
value v = 0.335±0.034. This suggests that there is only one critical state at the band center 
with all other states localized in the system on the infinite square lattice, in agreement with 
findings of Ref. |32|,20|. One parameter universal function A(x) is calculated and another 
discontinuity found, since A(x(E — > 0)) ^ A(E = 0, M). 

A puzzling feature of the critical exponent v of the ABD as well as of the RDF model 
fT7| , [T8|| is their apparent disagreement with the rigorous theorem of Chayes and Chayes et. 



al [43 , which states that v > 1 in two-dimensional quantum disordered systems in general. 



This question, however, requires further study and will be addressed elsewhere. 
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E E E 

FIG. 1. Density of states Pl{E) of the ABD model for w = 1 near the band center for system 
sizes L = 10, ...,60 in log-log plot. Full lines are fits to the Gade's form (|6|), while dotted lines are 
fits to the powerlaw. Fit to the Gade's form is done for all the points except the one closest to the 
band center, while fit to the powerlaw is done for all the bins except the three bins closest to the 
band center, which are the L-dependent parts of pi(E). 
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FIG. 2. Density of states of the ABD model for w = 1 near the band center. The graph is 
obtained from the data in Fig. |l|by removing the three bins closest to the band center, leaving the 
L-independent p(E). The fit is p(E) = Cexp (-KyJ- In \E\} /y/\E\, with k = 1.345 ± 0.005 and 
C = 1.30 ±0.03. 
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FIG. 3. Distributions of level spacings (after unfolding of the spectrum) Di(s) between i and 
i + 1 level, counted from the band center. D\ (s) is distinctly different than other distributions due 
to the presence of the symmetry. 
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FIG. 4. Dependence of the average IPN for several different q's and L = 80 on the bin size AE. 
The vertical dashed lines represent energies Ei(V)j1 (see text for discussion). 
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FIG. 5. Dependence of the average IPN for several different L and q = 0.9 on the bin size AE. 
The vertical dashed lines represent energies E2(L)/2 (see text for discussion). 
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FIG. 7. Scaling exponent of IPN, r q , and the goodness of fit coefficient 7 ?; calculated from the 
data in Fig. ^ (except L = 40 results). 
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FIG. 8. Multifractal spectrums D q at the band center and nearby energies. 
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FIG. 10. The largest renormalized inverse Lyapunov exponents A(E,M), calculated for long 
stripes of width M and various energies near the band center (full line). The energies range from 
log \E\ = — 1, — 1.5, — 5, and additional logE = —6. Dashed lines connect points of the second 
largest ILE, A 2 (E, M), for energies E = 10" 5 , 1(T 6 , and M < 20. 
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FIG. 11. One parameter universal function A(M/£(E)) for the ABD model. The four separate 
points are A(E = 0, M) for M = 16, 32, 64, 128, corresponding to the band center critical state. The 
inset shows the calculated localization length £(E) in units of the localization length £(E = 0.1), 
with error bars smaller than the symbol size. 
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